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Minkowski's theory of moving media is extended hereby to the anisotropic case. The corre- 
sponding Maxwell-Minkowski equations have been derived under the condition that the velocity of 
the moving medium is small compared to the velocity of light. As an application of that theory, it 
is shown that the characteristics of a plane wave propagating in a drifting magneto-ionic plasma can 
conveniently be determined from the constitutive parameters of the plasma without drifting. The 
equivalence between the convection current model adopted by Bell and Helliwell and the polarization 
current model suggested by Unz and correctly interpreted by Epstein, Bell, Smith, and Brice is also 
pointed out. 



1. Introduction 

The foundation of the electrodynamics of moving 
media was laid by Minkowski [1908] and based upon 
the special theory of relativity. For a moving iso- 
tropic medium, the complete treatment is given by 
Sommerfield [1952]. When the velocity of the 
moving isotropic medium is small compared to the 
velocity of light, the resultant wave equations can be 
greatly simplified [Tai, 1964b]. Several technical 
problems arising from these wave equations have 
already been investigated [Compton and Tai, 1964; 
Collier and Tai, 1964; and Tai, 1964a]. 

In this report we shall present an extension of Min- 
kowski's theory to moving anisotropic media. The 
theory is finally applied to a drifting magneto-ionic 
medium to determine the normal modes of plane 
waves which can exist in such a medium. The char- 
acteristic equation for the index of refraction so ob- 
tained by the first order relativistic transform method 
is the same as the one first derived by Bell and Helli- 
well [1959] using the convection current model and 
later by Epstein, Bell, Smith, and Brice [1963] using 
the polarization current model suggested by Unz 
[1962]. The equivalence between these models can 
be traced through Minkowski's theory. 



1.1. First-Order Relativistic Transformation of Field 
Vectors 

In this section, we shall first review the relativistic 
transformation of the field vectors defined in Max- 
well's equations as first recognized by Einstein and 
later elaborated by Minkowski. Two inertial systems 
(jc, y, z) and (x\ y' , z') are in relative motion. The 
primed system is assumed to be moving with a velocity 
v with respect to the unprimed system. The time vari- 
able in the two systems will be denoted, respectively, 
by t and t' . Maxwell's equations and the equation 
of continuity defined in the systems are 



V'X£": 
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If the velocity v is small compared to the velocity of 
light, the first-order relativistic transformation be- 
tween the two sets of field vectors, resulting from 
neglecting terms of the order of (v/c) 2 > is given by 
[Cullwick, 1957] 



D + - 9 vXH = e' -(E+vXB) 



B — -vXE = fjL f (H-vXD) 



(21) 
(22) 



^here 



E'=E+UxB 
B'=B-\vXE 



H'=H-vxD 



J'=J-pv 

l ~ 1 
P =p-^>v'J 



c 2 =lljJLoe . 



J-pv = o-' -(E+vXB) 



(23) 



(7) 

By eliminating D or B from (21) and (22) and neglect- 
(8) ing terms of the order of {v/c) 2 , we have 



(9) 
(10) 
(11) 

(12) 

(13) 



If two material field vectors P and M are used, such 
that one writes 



D = e<,E+P 



B = fjL (H + M) 



(14) 



(15) 



and similarly for the primed quantities, then the first- 
order relativistic transformation of the material field 
vectors can be deduced from (7-10). They are 



p'=p-l 



vXM 



M' = M + vXP. 



(16) 



(17) 



1.2. First-Order Constitution Relations Between the 
Field Vectors Defined in the Unprimed System 

Let us assume that constitutive relations for an 
anisotropic medium in the primed system are known. 
In the most general case, they are 



D' = 



E' 



B'=il'-H' 



(18) 



(19) 
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~E + vx(jl' 


H) 


-\vXH 


(24) 


B=ji'- 


H-vX(7 


■E) 


+ ^vXE 


(25) 


J=p- 


E + vx(jl' 


■H) 


+ pv. 


(26) 



Since V ■ D = p, (26) may be written in the form 



y=o-' 



E + v X (n' ■ H) 



+ »V •(€'•£). (27) 



Again, terms of the order of (v/c) 2 have been neglected 
in deriving (27). Equations (24), (25), and (27) contain 
the first-order constitutive relations between the field 
quantities defined in the unprimed system in terms 
of the constitutive parameters known or given in the 
primed system. By substituting these equations into 
(4) and (5), we obtain the Maxwell-Minkowski wave 
equations for the field vectors E and H in a moving 
anisotropic medium which represent a formal exten- 
sion of Minkowski's theory of moving media to the 
anistropic case. 

2. Magneto-Ionic Theory for Drifting Plasma 

In this section, we shall apply the above formulation 
to derive the characteristic equation for the index of 
refraction for a plane wave propagating in a drifting 
magneto-ionic medium. The convection current 
model will be treated first. In such a model the 
plasma is assumed to consist of free electrons with 
a drift velocity ¥. The primed system is attached to 
the drifting electrons and the unprimed system is 
fixed to an observer. According to the well-known 
theory of a stationary magneto-ionic jnedia, the rela- 
tion between the convection current J' and a harmon- 
ically oscillating, electric field E' with frequency a>' 
can be written in the form 



J' = <t'-E', 



(20) 



where the sign " = " is used to denote a dyadic. _ If the 
medium is dispersive, it is understood that "e', Ji f , and 
W would be functions of the frequency defined in the 
primed system, which will be denoted by (x)' '. By sub- 
stituting (7) to (11) into (18) to (20), we obtain 



>3'=E' 



where the dyadic r' is given by 



(28) 



X 



(o e 



(l-jZ')I+jY' 



(29) 
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The parameters contained in (29) are defined as follows: plane wave to be propagating in the z-direction in the 

unprimed system so that all the field components 



V 



,_ N\e[ 2 _,*>, 



Ffl 



(x) f , = plasma frequency 



Z' 



(O c 



(o' c = collision frequency, 



Y' = 





■y: 







r = Yg+Y$ + Y$-- 



Y'y 

-Y* 





fio\e\Ho 
mco' 



ft)' 



co„, = gyromagnetic frequency. 

The dyadic r' is clearly the recriprocal of tr' previously 
introduced in (20), i.e., 



/. 



In the convection current model, the medium is con- 
sidered to be unpolarized and unmagnetized, hence, 



6'= Co/ 

fJi' =fJbJ. 

Equations (24), (25), and (27), therefore, reduce to 

D = €of 

B = fJLoH 

7' • [J-vV • (iu&] = E + iiovxH. 



(31) 
(32) 

(33) 

(34) 
(35) 



On account of (33) and (34), (4) and (5) become, for a 
harmonically oscillating field, 



V X E =—j(DfjL H 
VxH=j+jcoeoE. 



(36) 
(37) 



Equations (35)_through (37) are the three basic equa- 
tions relating E, H, and J in a drifting magneto-ionic 
medium. They have been derived by applying 
Minkowski's first-order relativistic transform method. 
To derive the characteristic equation for the index 
of refraction, let us assume a harmonically oscillating 



have a dependence of e> {(1}t ~ kz) . The wave number k, 
the index of refraction /i, the phase velocity V p and 
the angular frequency (o are related as follows: 






(38) 



As a result of the assumption of a plane wave propa- 
gating in the z-direction, (36) and (37) can be written as 

-jkz XE = -jcofioH (39) 

-jkz XH = J +ja)€oE. (40) 

Eliminating H between (39) and (40), we obtain 

J = -jwe,[(l - n 2 )E + n*E z z\. (41) 

Substituting the expressions for H and J as given by 
(39) and (41) into (35), we obtain a homogeneous equa- 
tion for E 

-jcoP • [( 1 - n 2 )E + (n 2 2 - np)E z ] = E + npX(zXE) 

(42) 



(30) where 



P = v/c. 



(43) 



To find the characteristic equation for n, the primed 
quantities contained in~r' must first be converted into 
explicit functions of n. Because of the Doppler shift, 
the relationship between co' and co, accurate to the 
order of v/c, is given by 



= w(l- 



■V Z ) 



-(x)(l V Z ) 

C 



'Ml-nPzh 



(44) 



The parameters contained in r', therefore, can be 
written as: 



W) (l-np z f\a> 



\ 2 = * 
a)' l-np z \o>) 1-nfig 



r 



0)' 



l-np z 



1-nPz 



(45) 
(46) 

(47) 
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hence, 



. ,=, l-nfe 



-jwe r 



X 



1-7 



* ^7+y ' 



l-n# 



l-«& 



=j(l-n^-;Z)/+ 7 T. 



(48) 



On account of (48), the explicit expression for (42) is 
given by 

[Q.-np z -jZ)T+j?\ ■ [(l-n?)E + (n 2 z-n]})Ez] 

=X[E + nJBX(zXE)]. (49) 

Equation (49) can be written in the form 

2 1=0. (50) 

Where the element of the dyadic A are found to be 
A xx = (n 2 - 1) (1 - npz -JZ) + (l-n/3 z )X 
A xy =-j(n*-l)Y z 

A xz = (1 - nfi z -jZ)np x -jnpyYz -j(l - np z )Y y 
A yx =j(n*-\)Y Z 

A yy = (rc 2 - 1) (1 - np z -jZ) + (1 - np z )X 
A vz = (1 - nB* -yZ) RJ 8, +jnp x Y z +jd ~ nB z )Y x 
A zx =-j(n 2 -l)Y y + nB x X 
A zy =j(n i -l)Y x + nByX 

A zz =-(l-nB z )(l- n/3 z -jZ) -jnfrrYy +jnB y Y x + X. 

(51) 

By setting the determinant of A equal to zero, one 
obtains the characteristic equation for n. Equation 
(50) is the same as the one first obtained by Bell and 
Helliwell [1959] 2 if we discard in their work the part 
contributed by the ambient plasma. It is also identical 
to the one obtained by Epstein, Bell, Smith, and Brice 
[1963] after correctly interpreting the work of Unz 
[1962]. The alternative derivation presented here 
appears to provide a clearer view as to how the original 
Appleton-Hartree matrix was transformed as a result 
of the motion. It also shows another application of 
Minkowski's powerful theory. While the treatment of 
Bell and Helliwell may be considered as "classical" 
if one accepts the Lorentz force equation as a prere- 



2 Bell and Helliwell considered only the case where v and H have the same direction. 
This restriction can easily be removed. They also included a stationary plasma in addition 
to the drifting one. The part contributed by the stationary plasma can be included in our 
treatment if necessary. There are some minor algebraic mistakes in Bell and Helliwell's 
work. For example the matrix for J MU equation (12), contains a parameter p t which is 
defined as klk\Y L that appears to be incorrect. The parameter / 2 defined in the list of sym- 
bols does not seem to fit the dimensional analysis. 



lativistic postulate, as we normally do, instead of as a 
consequence of the relativistic transformation, the 
treatment of Unz, Epstein, et al., is not. This is be- 
cause the E and H fields in the EPHMv formulation 
are only partial fields [1964b]. The force equation 
in that formulation therefore does not have the same 
structure as the Lorentz force equation. In order to 
justify that treatment, we shall treat the same problem 
based upon the polarization current model and point 
out its implications. 

3. Polarization Current Model 

In the book by Ratcliffe [1959], the Appleton- 
Hartree equation for a stationary magneto-ionic me- 
dium was derived by considering the medium to be 
made of polarized matter or bound electrons. We 
shall call such a model "the polarization current 
model" in contrast to the convection current model 
presented in the previous section. In the polariza- 
tion current model, one has 



/'=0 



M' = 



(52) 



(53) 



and 



P' = eox'E' (54) 

where, numerically, the dyadic susceptibility \' ^ s 
related to <x' defined by (30) in the following way 



j(o eox = <r 



As a result of using P\ (24) and (25) become 
D=7' '{E-^fjio^xH)-\vxH 



B = /JL 



H-vX(e' -E) 



+ \vXE 



where 



e' = e (/ + x')- 
Two alternative expressions for (56) and (57) are 

Z> = 6o£; + €oX' ■(E + IXoVXlf) 



B = (j-oH — ijlov X (eox' • E). 



If we denote 



eox'-(E + HovXH)=P 



(55) 

(56) 
(57) 

(58) 

(59) 
(60) 

(61) 
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then 



If we let 



B = tio(H-vXP). 



(62) 



Within the order of accuracy of v/c, the apparent dif- 
ference between (60) and (63) is of no consequence. 
The basic equations involved in the polarization current 
model are summarized below: 



P = -NeR n 



(75) 



(63) then, (74) is equivalent to 



d 2 R„ ,dR p 



nD _ 

-—e(E + n<r£? X H t) + ^v X H p ). 



(76) 



V X E = -ja>no(H -vXP) 

V XH =jo)(eoE + P) 

P = e x' -(E + ^vXH). 



(64) 
(65) 
(66) 



These equations play a role similar to (35) through 
(37) in the convection current model. As far as the 
result is concerned it will be shown that two models 
do yield the same answer as far as the characteristic 
equation for n is concerned. Equations (64) through 
(66) may be transformed into another form by intro- 
ducing a field vector H p defined by 



H f) = H-vXP 



that gives 



V XE^—jlO/JLoHp 

V x H p *= joKcoE + P) - V x ( v x P) 

P = eox' '(E + VsovxHpY- 



rn 



(68) 



(69) 



(70) 



A term of the order of (v/c) 2 is discarded in converting 
(66) into (70). Equations (68) and (69) are identical in 
form to the two equations used by Unz [1962] who 
originally investigated that problem based upon the 
so-called EHPMv formulation. Equation (70) can 
be transformed to the force equation interpreted by 
Epstein [1963]. Thus, if we write (70) in the form 



or 



J0)'P = (T' '(E+fJLoVXHp) 



j<o'r'-P = (E+tiovXH p ). 



(71) 



(72) 



Using r' given by (29) and substituting it into (72) we 
may reconstruct the differential equation correspond- 
ing to (72). We consider 



£[<!-*' 



which is the same as 



)I+jT 



P = E + fi.<,vXH p (73) 



OTO» 



N\ ( 



('-£) 



a) I ma) 



■E+/jlqvXH p . 

(74) 



Since t' is defined in the primed system, the partial 
derivation of a function with respect to t' is the same 

as the material derivative — = \-(v • V) evaluated in 

at at 

the unprimed system as a result of the first-order 
Lorentz transformation, i.e., neglecting terms of the 
order of (v/c) 2 . Equation (76), therefore, is identical 
to the force equation in the EHPMv formulation. As 
we have mentioned earlier, as far as the character- 
istic equation for n is concerned, both the convection- 
current model and the polarization model provide the 
same answer. However, it is obvious from this dis- 
cussion that there are several delicate concepts in- 
volved in the polarization current model. It appears 
that these concepts and their acceptance can best be 
understood and justified with the aid of Minkowski's 
theory. Finally, it should be pointed out that although 
the classical method provides the correct answer 
when {v/c) 2 << 1 it fails when v is comparable to c. In 
that case, one must depend upon Minkowski's exact 
theory to formulate such a problem. 
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